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Abstract 



Numerically we solve the microscopic deterministic equations of motion with 
random initial states for the two-dimensional cfi 4 theory. Scaling behavior of 
the persistence probability at criticality is systematically investigated and the 
persistence exponent is estimated. 
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Recently the persistence exponent has attracted much attention. This exponent was first 
introduced in the context of the non-equilibrium coarsening dynamics at zero temperature 
P|||. It characterizes the power law decay of the persistence probability that a local order 
parameter keeps its sign during a time t after a quench from a very high temperature to 
zero temperature. For critical dynamics, the local order parameter (usually, a spin) flips 
rapidly and the persistence probability does not obey a power law. In this case, however, 
the persistence exponent 9 P can be defined by the power law decay of the global persistence 
probability p(t) that the global order parameter has not changed the sign in a time t after 
the quench from a very high temperature to the critical temperature ||, 

p(t) ~ & . (1) 

An interesting property of the persistence exponent is that its value is highly non-trivial 
even for simple systems. For the quench to zero temperature, for example, 9 P is apparently 
not a simple fraction for the simple diffusion equation and the Potts model in one dimension 
For the quench to the critical temperature, it is shown that the persistence exponent 
is generally a new independent critical exponent, i.e. it can not be expressed by the known 
static exponents, the dynamic exponent z and the recently discovered exponent 9 [|3]||. This 
relies on the fact that the time evolution of the global magnetization is not a Markovian 
process. Recent Monte Carlo simulations for the Ising and Potts model at criticality support 
the power law decay of the global persistence probability and detect also the non-Markovian 
effect J70. 

Up to now, the persistence exponent has been studied only in stochastic dynamic systems, 
described typically by Langevin equations or Monte Carlo algorithms. From fundamental 
view points, both equilibrium and non-equilibrium properties of statistical systems can be 
described by the microscopic deterministic equations of motion (e.g. Newton, Hamiltonian 
and Heisenberg equations) even though a general proof does not exist. With recent de- 
velopment of computers, gradually it becomes possible to solve numerically the microscopic 
deterministic equations of motion. For example, the O(N) vector model and the XY model 
have been investigated pHIIf. The results confirm that the deterministic equations describe 
correctly second order phase transitions. The static critical exponents are estimated and 
agree with existing values. More interestingly, in a recent paper short-time dynamic behav- 



ior of the deterministic dynamics starting from random initial states has been studied |L2 
The short-time dynamic scaling was found and the estimated value of the dynamic exponent 
z is the same as that of the Monte Carlo dynamics of the Ising model. 

The purpose of this letter is to study the critical scaling behavior of the global persistence 
probability and measure the persistence exponent in microscopic deterministic dynamic 
systems, taking the two-dimensional 4 theory as an example. 

The Hamiltonian of the two-dimensional 4 theory on a square lattice is written as 



(2) 



Langevin equations at zero temperature are also deterministic, but they are at mesoscopic level 
and generally different from the microscopic deterministic equations of motion. 
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with 7Tj = <pi and it leads to the equations of motion 



1 

& = + <l>i-n ~ 2( t>i) + m2< Pi ~ ■ (3) 

Energy is conserved during the dynamic evolution governed by Eq. ([J). As discussed in 
Refs. p,12|, a microcanonical ensemble is assumed to be generated by the solutions. In this 



case, the temperature could not be introduced externally as in a canonical ensemble, but 
could only be defined internally as the averaged kinetic energy. In the dynamic approach, 
the total energy is actually an even more convenient controlling parameter of the system, 
since it is conserved and can be input from the initial state. 

From the view point of ergodicity, to achieve a correct equilibrium state the microscopic 
deterministic dynamic system should start from a random initial state. Interestingly, this 
is just similar to the dynamic relaxation in stochastic dynamics after a quench from a very 
high temperature. Therefore, similar dynamic behavior may be expected for both dynamic 
systems. 

The order parameter of the 4 theory is the magnetization M(t) = Y^,i4>i{t)/L 2 with L 
being the lattice size. In this paper, we are interested in the global persistence probability 
p(t) at the critical point, which is defined as the probability that the not averaged order 
parameter has not changed the sign in a time t starting from a random state with small 
initial magnetization m . 

Following Ref. P,|T2||, we take parameters m 2 = 2. and A = 0.6 and prepare the initial 



configurations as follows. For simplicity, we set initial kinetic energy to be zero, i.e. (pi(0) = 
0. We fix the magnitude of the initial field to be a constant c, |0i(O)| = c, and then randomly 
give the sign to <pi(0) with the restriction of a fixed magnetization in unit of c, and finally 
the constant c is determined by the given energy. 

To solve the equations of motion @ numerically, we simply discretize (pi by ((pi(t + At) + 



4>i(t - At) - 2(pi{t)) / (At) 2 . According to the experience in Ref. ||12|| , At is taken to be 0.05. 
After an initial configuration is prepared, we update the equations of motion until the mag- 
netization changes its sign. The maximum observing time is t = 1000. Then we repeat the 
procedure with other initial configurations and measure the persistence probability p{t). In 
our calculations, we use fairly large lattices L = 128 and 256 and samples of initial configu- 
rations range from 3 500 to 30 000 depending on initial magnetization mo and lattice sizes. 
The smaller mo and lattice size L are, the more samples of initial configurations we have. 
Errors are simply estimated by dividing total samples into two or three subsamples. Com- 
pared with Monte Carlo simulations, our calculations here are much more time consuming 
due to the small At. 

According to analytical analyses and Monte Carlo simulations in stochastic dynamic 
systems, at the critical point and in the limit mo = 0, p(t) should decay by a power law as 
in Eq. (|1]). Our first effort is to investigate whether in microscopic deterministic dynamics 
p(t) obeys also the power law and measure the persistence exponent 8 p . 

Here we adopt the critical energy density e c = 21.1 from the literature [§,0- In Fig. [IJ, 
the persistence probability p(t) is displayed on a log-log scale for lattice sizes L = 256 and 128 
with solid lines and a dash line respectively. For L = 256 simulations have been performed 
with two values of initial magnetization mo = 0.003 and 0.0015. These straight lines convince 
us the power law behavior of p{t). Skipping data within a microscopic time scale t ~ 100, 
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from the slopes of the curves one estimates the persistence exponent 9 P = 0.252(6) for both 
values of m . This shows that there is already no effect of finite m . The curve for L = 128 
and m = 0.0015 is roughly parallel to that of L = 256. One measures the the slope 
9 P = 0.251(1) in the time interval [100, 500] but 9 P = 0.232(10) in [100, 1000]. This indicates 
that some finite size effect exists still for L = 128 after t = 500 but is negligible small for 
L = 256. 

If the time evolution of the magnetization is a Markovian process, from theoretical view 
points the persistence exponent will be not an independent critical exponent and it will take 
the value a p , which relates to other exponents through 

a p = -9 + (d/2 - P/u)/z . (4) 

In Table [| values of the exponent 8 P , z, 9 and a p for the 4 theory are given in comparison 
with those of the kinetic Ising model induced by local Monte Carlo algorithms. As is the 
case of the Ising model, the exponents 9 P and a p for the <j) A theory differ also by about 10 
percent. This represents a rather visible non-Markovian effect in the time evolution of the 
magnetization. 

For equilibrium states, it is generally believed that the 4 theory and the Ising model are 
in a same universality class. Results from numerical solutions of the deterministic equations 
also support this ||. From the short-time dynamic approach [|12[], within statistical errors 
the dynamic exponent z for the microscopic deterministic dynamics of the <f) 4 theory is the 
same as that of the kinetic Ising model with Monte Carlo algorithms but the exponent 9 
differs by several percent. In Table |, we see that 9 P for the </> 4 theory is also several percent 
bigger than that of the Ising model. However, by feeling we still think that the 4 theory 
and the Ising model are very probably in a same persistence universality class. These some 
percent differences of the exponents come probably from that the critical point e c has not 
been very accurate or there are some corrections to scaling f] and uncontrolled systematic 
errors. Actually, we will see below that the critical energy density could be somewhat lower 
than e c = 21.1 and it would yield a slightly smaller 9 P . 

Our second step is to investigate the scaling behavior of the persistence probability in 
the neighborhood of the critical energy density. From general view points of physics, one 
may expect a following scaling form 

p(t,r) =r ep F(t 1/uz T) . (5) 

Here r = (e — e c )/e c is the reduced energy density. When r = 0, the power law behavior in 
Eq. ([TJ) is recovered. When r differs from zero, the power law will be modified by the scaling 
function F(t 1 ' l/Z T). In principle, this fact may be used for the determination of the critical 
energy density. In Fig. |2|, p(t, r) for L = 256 and m = 0.003 is plotted for three different 
energy density e = 20.9, 21.1 and 21.3. In the figure, we see that the solid line shows the 



2 In deterministic dynamics, energy is conserved and it couples to the order parameter. Therefore, 
it is believed that the deterministic dynamics belongs to the dynamics of model C rather than 
model A. Standard local Monte Carlo dynamics of the Ising model is dynamics of model A. In two 
dimensions, model A and C are the same but maybe up to a logarithmic correction [13,1^]. 
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best power law behavior among the three curves. However, an very accurate estimate of e c 
could not be achieved so easily, since p(t, r) is not so sensitive to the energy. According to 
our data, we estimate e c = 21.06(12). Within errors, it is consistent with e c = 21.1 given in 



@ and e c = 21.11(3) in fl2]| . We should point out, that the exponent 9 P will be 0.245, if it 
is measured at e c = 21.06. It is closer to that of the kinetic Ising model, as discussed above. 

In order to have more understanding of the scaling form @, we differentiate with respect 
to the energy density on both sides of the equation and obtain 

d T lnp{t,r) | T=0 ~t lA " • (6) 

Using the data of Fig. 0, we can approximately calculate d T \np(t, t)| t=0 and the result is 
displayed in Fig. || Even though there are some fluctuations, power law behavior is still 
seen. The best fitted slope of the curve gives 1/uz = 0.47(4) in the time interval [100, 1000]. 
Taking z = 2.15(2) as input, one obtains v = 0.99(8). Compared with v — 1 for the Ising 
model, this result supports that the 4 theory with deterministic dynamics and the Ising 
model are in a same universality class. 

Finally we study the scaling behavior of the persistence probability in case the initial 
magnetization is not so small and its effect can not be neglected. Following Ref. ||, we 
assume a finite size scaling form 

p(t, L, m ) = r dp F(t 1/z L~\ t XOp/z m ) . (7) 

Here the energy density has been set to its critical value and xq p is the scaling dimension 
of the initial magnetization mo- It was discovered that for the Ising model with Monte 
Carlo dynamics, the value x$ p = 1.01(1) is 'anomalous', i.e. it is different from the scaling 
dimension of the initial magnetization x = 0.536(2) measured from the time evolution of 
the magnetization or auto-correlation ||. The origin should be that p(t, L, mo) is a non-local 
observable in time t. It remembers the history of the time evolution. 

To verify the scaling form (0) and estimate xo p , we perform a simulation with the lattice 
size Li = 256 and initial magnetization m 01 = 0.0151. Suppose the scaling form (j^j holds, 
one can find an initial magnetization m 02 with the lattice size L 2 = 128 such that the curves 
of p(t, L, mo) for both lattice sizes collapse. Practically we haved performed simulations for 
L 2 = 128 with two initial magnetizations, mo = 0.0272 and 0.0349. By linear extrapolation, 
we obtain data with m between 0.0272 and 0.0349. Searching for a curve best fitted to 
the curve for Li = 256, we determine mo2- In Fig. |], such a scaling plot is displayed. The 
lower and upper solid lines are the persistence probability for L 2 = 128 with mo = 0.0272 
and 0.0349 respectively, while the dashed line is the properly rescaled one for L\ = 256 with 
moi = 0.0151. The solid line fitted to the dashed line represents the persistence probability 
for L 2 = 128 with mo2 = 0.0313(3). Since the microscopic time scale is t m i C ~ 100, nice 
collapse of the two curves can be observed only after t ~ 100. From the scaling form (^), 
m 2 = 2 XOp m i and one estimates x 0p = 1.05(1). This value is very close to x 0p = 1.01(1) 
for the kinetic Ising model. 

In conclusions, we have numerically solved the microscopic deterministic equations of 
motion with random initial states for the two-dimensional <f) A theory and systematically 
investigated the critical scaling behavior of the persistence probability. As summarized in 
Table |, the estimated exponents 6 P and x 0p are very close to those of the kinetic Ising model 
induced by local Monte Carlo algorithms. What would be the dynamic and static behavior 
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of the deterministic dynamics starting from more general initial states is an interesting work 
in future. 

Acknowledgements: Work supported in part by the Deutsche Forschungsgemeinschaft 
under the project TR 300/3-1. 



6 



REFERENCES 



[1] B. Derrida, A.J. Bray and C. Godreche, J. Phys. A27, L357 (1994). 

[2] A. J. Bray B. Derrida and C. Godreche, Europhys. Lett. 27, 175 (1994). 

[3] S.N. Majumdar, A.J. Bray, S. Cornell and C. Sire, Phys. Rev. Lett. 77, 3704 (1996). 

[4] B. Derrida, V. Hakim and V. Pasquier, Phys. Rev. Lett. 75, 751 (1995). 

[5] S.N. Majumdar, C. Sire, A.J. Bray and S. Cornell, Phys. Rev. Lett. 77, 2867 (1996). 

[6] K. Oerding, S.J. Cornell and A.J. Bray, Phys. Rev. E56, R25 (1997). 

[7] D. Stauffer, Int. J. Mod. Phys. C7, 753 (1996). 

[8] L. Schiilke and B. Zheng, Phys. Lett. A 233, 93 (1997). 

[9] L. Caiani, L. Casetti and M. Pettini, J. Phys. A31, 3357 (1998). 
[10] L. Caiani, L. Casetti, C. Clementi, G. Pettini, M. Pettini and R. Gatto, Phys. Rev. 

E57, 3886 (1998). 
[11] X. Leoncini and A.D. Verga, Phys. Rev. E57, 6377 (1998). 
[12] B. Zheng, M. Schulz and S. Trimper, Phy. Rev. Lett. 82, 1891 (1999). 
[13] PC. Hohenberg and B.I. Halperin, Rev. Mod. Phys. 49, 435 (1977). 
[14] K. Oerding, private communication, 1999. 
[15] B. Zheng, Int. J. Mod. Phys. B12, 1419 (1998), review article. 



7 
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FIG. 1. The persistence probability is displayed in log-log scale for e c = 21.1. Solid lines and a 
dashed line are for lattice sizes L = 256 and 128 respectively. 
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FIG. 2. The persistence probability is displayed in log-log scale for the energy densities e = 20.9, 
21.1 and 21.3 (from above). 
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FIG. 3. The logarithmic derivative of the persistence probability with respect to the energy 
density is plotted in log-log scale, using the data in Fig. |^. The dashed line is the best fitted 
straight line 
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FIG. 4. The scaling plot for the persistence probability. The lower and upper solid lines are the 
persistence probability for Li = 128 with tuq = 0.0272 and 0.0349 respectively, while the dashed 
line is the properly rescaled one for L\ = 256 with moi = 0.0151. The solid line fitted to the 
dashed line represents the persistence probability for L2 = 128 with mo2 = 0.0313(3). 
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TABLES 





e p 


z 


9 


a p 


1/uz 


V 


XOp 




0.252(6) 


2.15(2) 


0.176(7) 


0.231(7) 


0.47(4) 


0.99(8) 


1.05(1) 


Ising 


0.238(3) 


2.155(3) 


0.191(1) 


0.215(1) 




1 


1.01(1) 



TABLE I. The critical exponents measured for the </> 4 theory in comparison with those of the 
Ising model. The value of v for the Ising model is exact, while others are taken from Table 2 in 
Ref . [|l5[] and from Ref . ||] . z and 9 for the 4 theory are from Ref . |li| . To calculate a p for the </> 4 
theory /3/u = 1/8 is taken as input. 
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